Nuclear modification of high-p T hadron spectra in p + A collisions at LHC 
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Multiple parton scatterings in high-energy p + A collisions involve multi-parton correlation inside 
the projectile and color coherence of multiple jets which will lead to nuclear modification of final 
hadron spectra relative to that in p + p collisions. Such modifications of final hadron spectra in 
p + A collisions are studied within the HIJING 2.1 model. Besides parton shadowing and transverse 
momentum broadening, which influence final hadron spectra at intermediate pr, modification of 
parton flavor composition due to multiple scattering and hadronization of many-parton jets system 
are shown to lead to suppression of final hadrons at large pr in p + A collisions at the Large Hadron 
Collider. 

PACS numbers: 24.85.+p,25.75.-q,25.75.Bh 



In search for the formation of quark-gluon plasma 
(QGP) in high-energy heavy-ion collisions, jet quenching 
or suppression of high-pT hadrons [l[ can be used for to- 
mographic study of medium properties [!, H|. However, 
initial multiple parton scatterings in cold nuclei before 
the formation of QGP can also lead to nuclear modifi- 
cation of final high-pT hadron spectra that needs to be 
understood for a quantitative study of QGP properties 
through jet quenching. 

Parton shadowing or nuclear modification of parton 
distributions inside large nuclei is one of the most stud- 
ied cold nuclear effects. The QCD evolution, however, 
dramatically reduces the effect on final high-pf hadron 
spectra in p + A collisions [iHZj]. Modification of par- 
ton distributions in the projectile nucleon through par- 
ton energy loss || |j| is estimated to have a small ef- 
fect on hadron spectra in the central rapidity region of 
high-energy p + A collisions. These cold nuclear effects 
have also been studied within the Color Glass Conden- 
sate (CGC) model fll)14l2^ . However, multiple parton 
correlations inside the projectile nucleon should also be 
considered. Energy-momentum and valence quark num- 
ber conservation alone could modify the momentum and 
flavor dependence of final-state parton and hadron spec- 
tra at large px- Hadronization of multiple-jet systems 
from multiple scattering in p + A collisions could also 
affect final hadron spectra at intermediate and large pr- 

In this paper, we study the nuclear modification of 
parton and hadron spectra due to parton shadowing, 
initial and final-state transverse momentum broadening, 
final-state parton flavor composition and hadronization 
of multiple jets within the HIJING [H| Monte Carlo 
model. It is based on a two-component model for hadron 
production in high-energy nucleon and nuclear collisions. 
The soft and hard components are separated by a cut-off 
Po in the transverse momentum exchange. Hard parton 
scatterings with px > Po are assumed to be described 
by the perturbative QCD (pQCD) and soft interactions 
are approximated by string excitations with an effective 



cross section u so f t . In HIJING 2.0 [14(, parton distri- 
bution functions (PDF's) are given by the Gluck-Rcya- 
Vogt(GRV) parameterization [15(. An energy-dependent 
cut-off po(v^) and s °ft cross section f7 so f t (-y/s) are as- 
sumed in order to fit experimental values of the total 
and inelastic cross sections of p + p/(p) collisions and 
the energy-dependence of the central rapidity density of 
charged hadron multiplicities. HIJING 2.0 can describe 
most features of hadron production in p + p collisions 
at colliding energies up to 7 TeV at the Large Hadron 
Collider (LHC) Q. 

In p + A and A + A collisions, one has to consider 
nuclear modification of initial parton distributions due to 
multiple scatterings. HIJING 2.0 employes a factorized 
form of parton distributions in nuclei, 

f a /A(x, Q 2 ,b) = AR^(x, Q 2 ,b)f a/N (x, Q 2 ), (1) 

where Q is the momentum scale, b the impact-parameter, 
x the momentum fraction carried by the parton a and 
Ra(x,Q 2 ,b) is the impact-parameter-dependent nuclear 
modification factor as given by the new HIJING parame- 
terization [l6|. It was fit to the centrality dependence of 
the pseudo-rapidity multiplicity density per participant 
pair 2dN c h/drj/N paI t at the Rclativistic Heavy-ion Col- 
lider (RHIC), and also describes the data at the LHC 
energies as well [13] . 

In addition to modification of initial parton distri- 
butions, multiple scatterings inside a nucleus should 
also lead to transverse momentum (far) broadening of 
both initial and final-state partons The fc<r-broadcning 
is responsible for the Cronin effect or enhancement of 
intermediate-pT hadron spectra in p + A collisions fl7j . 
It is implemented in HIJING 2.1 as a fey-kick by each 
binary nucleon-nucleon scattering to the initial partons 
of a hard scattering as well as final-state partons pro- 
duced prior to the current scattering. The fc^-kick for 
each scattering follows a Gaussian distribution and an 
energy dependence of the width, 
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FIG. 1: (color online) Charged hadron rapidity distribution 
from HIJING 2.1 with (solid) and without (dot-dashed) par- 
ton shadowing in p(d) + A and p + p (dashed) collisions at 
y/s NN = 200 GeV (upper) and 5.0 TeV (lower) as compared 
to data from the STAR pj| and PHOBOS experiment. 



is assumed to describe the measured Cronin effect in p+A 
collisions. After the fey-kick, longitudinal momenta of 
partons are reshuffled pair-wise between projectile and 
target partons in order to ensure four-momentum con- 
servation. Consequently, introduction of the fey-kick will 
also influence hadron rapidity distributions. We have to 
fine-tune the gluon shadowing parameter s g [l4l , [l6j again 
to describe experimental data on rapidity distributions. 

Shown in Fig.[T]are HIJING 2.1 results on rapidity dis- 
tributions of charged hadrons in minimum-biased d + Au 
(upper panel) and p+Pb (lower panel) collisions (solid) at 
\/s NN = 200 GeV and 5.0 TeV, respectively, as compared 
to d + Au data from the STAR experiment [18| . Here the 
gluon shadowing parameter is set to s g = 0.28. Also 
shown are HIJING2.1 results for p + p collisions (dashed) 
as compared to the PHOBOS data [lj| and p{d) + A col- 
lisions without parton shadowing (dot-dashed). Gluon 
shadowing in HIJING 2.1 is seen to suppress the hadron 
yield in mid and forward rapidity region, especially at 
the LHC energy. Measurement of the hadron rapidity 
distribution at LHC therefore would provide further con- 
straints on the gluon shadowing. 

The hard processes in HIJING are based on lowest or- 
der pQCD in which the single jet inclusive cross section, 
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FIG. 2: (color online) Nuclear modification factor for parton 
Pt spectra in p(d) + A collisions at \/s NN = 200 GeV (upper) 
and 5.0 TeV (lower) from HIJING 2.1 with (solid), without 
(dashed) parton shadowing and with de-coherent hard scat- 
terings (DHC) (dot-dashed and doted). 



/a/A (^2, Q 2 , b). Here xi )2 = 2?r(e ±yi +e ±V2 )/y/s are frac- 
tional momenta of the initial partons, y± t 2 are rapidities 
of the final-state parton jets and K ss 2 accounts for 
higher-order corrections. The nuclear thickness function 
is normalized to J d 2 btA(b) = 1. In the central rapid- 
ity region (2/1,2 = 0), transverse momenta of final-state 
jets determine the momentum fractions 2:1,2 = 2pr/\/s 
of the initial colliding partons. Therefore, px spectra 
of produced partons via hard scatterings in p + A colli- 
sions should depend on the nuclear modification of par- 
ton distributions in nuclei. To study parton shadowing 
and other cold nuclear effects, we calculate the nuclear 
modification factor for final-state parton and hadron px 
spectra, 



R p a(pt) 



dN pA /dyd 2 p T 
(N hin )dN pp /dyd 2 p T ' 



(4) 



is directly proportional to nuclear parton distributions 



where (iVbin) is the average number of binary nucleon- 
nuclcon interactions in p + A collisions. 

Shown in Fig. [2] are parton modification factors from 
the HIJING 2.1 for minimum-biased d + Au and p + Pb 
collisions at \/s NN = 200 GeV and 5.0 TeV, respec- 
tively. The enhancement of parton spectra at interme- 
diate pt is the signature effect of fey-broadening, which 
should disappear at large pr- In the HIJING model, 
hard scatterings are simulated each time the projectile 
parton interact with a target parton alongside soft inter- 
actions with the remnants of colliding nucleons. HIJING 
ensures energy-momentum conservation by subtracting 
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FIG. 3: (color online) Quark to gluon ratio as a function of 
Pt in p + p (dotted) p(d) + A collisions at the RHIC (upper) 
and LHC (lower) energies from HIJING 2.1 without (solid) 
and with (dashed) de-coherent hard scatterings (DHC). 



energy-momentum transfers in the previous hard and soft 
interactions from the projectile and thus restricting en- 
ergy available for subsequent hard collisions. Since the 
time scale of hard scatterings is much shorter than soft 
interactions, such a coupling between hard and soft in- 
teractions might not be physical. As an option in our 
study, we turn off such a coupling and do not restrict 
energy available for hard scatterings. We denote the re- 
sults as DHC (de-coherent hard scattering). One can 
see from Fig. [21 both parton shadowing and soft-hard 
coupling suppresses high-py parton production. These 
features of parton spectra will be reflected in the final 
hadron spectra after hadronization. 

In p + A collisions, the projectile proton scatters with 
multiple nucleons inside the target, each with finite prob- 
ability of hard parton scattering. The leading contribu- 
tion comes from independent hard scatterings, each in- 
volving independent initial partons from projectile and 
target nucleons. In principle, a fast projectile parton can 
also suffer multiple hard scatterings which will affect the 
final-state parton spectra but do not increase the pro- 
jectile's contribution to the final-state parton multiplic- 
ity. Momentum correlation of multiple partons inside the 
projectile is neglected beyond the conservation of the to- 
tal momentum. Flavor conservation, however, will limit 
the availability of valence quarks from the projectile for 
each of these multiple independent hard scatterings as 
the projectile plunges through the target nucleus. Such 
limitation on the number of valence quarks available for 
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FIG. 4: (color online) Proton to pion ratio at RHIC and 
LHC from HIJING 2.1 with different options: default (solid), 
without shadowing (ashed), DHC (dot-dashed), DHC with- 
out parton shadowing (dotted) in d + Au collisions at RHIC 
and p + Pb at LHC. The ratio in p + p collisions is also shown 
(long dashed). The RHIC experimental data are from STAR 



each independent hard scattering will lead to a change 
in the flavor composition of produced partons per aver- 
age binary nucleon-nucleon collision. Shown in Fig. [3] 
are the ratios of produced quark and gluon pt spectra 
in the central rapidity region of d + Au and p + Pb col- 
lisions at \/~s NN = 200 GeV and 5.0 TeV, respectively, 
as compared to that in p + p collisions. The quark num- 
ber conservation in multiple scatterings is seen to lead to 
smaller quark/gluon ratio in p + A relative to p+p colli- 
sions. The suppressed quark yield and high quark/gluon 
ratio, especially at RHIC, will lead to the modification of 
final-state parton spectra at large pt as shown in Fig. [2] 
Because gluon fragmentation functions are softer than 
quarks, the increased fraction of gluon jets in p + A colli- 
sions will lead to apparent suppression of high px hadron 
spectra relative to p + p collisions. 

The modified flavor composition of the final parton 
due to multiple scattering in pA collisions could also 
change the flavor composition of final hadrons due to 
flavor-dependence of the fragmentation functions. Gluon 
jets in general have higher proton-to-pion ratio than in 
quark jets. HIJING has a strong gluon shadowing which 
suppresses larger pt gluon jets. This should reduce the 
proton-to-pion ratio in HIJING 2.1 results as compared 
to the results without gluon shadowing as shown in Fig. [4] 
This reduction of proton-to-pion ratio by gluon shadow- 
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ing is much stronger at the LHC energy. Nevertheless, 
the proton-to-pion ratio in d + Au collisions at RHIC and 
p + Pb collisions at LHC will still be larger than that in 
p + p collisions because of the enhanced gluon-to-quark 
ratio in p + A collisions relative to that in p + p as shown 
in Fig. [3] Such a trend is already seen in the STAR data 
[2lj at RHIC. The flavor composition of final hadrons in 
p + A collisions could also be influenced by other non- 
perturbative mechanisms such as parton recombination 
during hadronization (22| and hadronic final-state rescat- 
tering [23| . 

Another process that can also modify final hadron 
spectra in p+A collisions is jet hadronization or fragmen- 
tation. HIJING follows the same routine as in PYTHIA 
that groups jet partons together with project and tar- 
get remnants to form color-singlet string systems. In 
HIJING, jet parton showers from all hard scatterings, 
including initial and final-state radiations, are ordered in 
rapidity and gluons are then connected to the valence 
quark and diquark of the associated projectile or target 
nucleons as kinks to form string systems. These string 
systems are hadronized according to the Lund string 
model [2(j. The projectile can undergo multiple scat- 
terings in p + A collisions. Its string system therefore 
has more gluons attached as compared to that in p + p 
collisions. Hadrons from such a string system are ex- 
pected to be softer than independent fragmentation of 
individual gluons. To test the sensitivity of nuclear mod- 
ification of final hadron spectra to jet hadronization, we 
also consider the independent fragmentation mechanism 
in HIJING for both p + p and p + A collisions. 

Shown in Fig. [5] are nuclear modification factors for 
charged hadrons in d + Au and p + Pb collisions at 
\/s NN = 200 GeV and 5.0 TeV, respectively, as com- 
pared to existing data from both the PHENIX and STAR 
cxpok, eriment at RHIC [U, [25|. At low p?, hadron 
production is dominated by hadronization of soft strings 
which is proportional to the number of participant nu- 
cleons Ap ar t = 1 + Abin- The nuclear modification factor 
should approach R p a ~ (1 + l/Abi n )/2 « 1/2 in the 
limit Abin 3> 1 • At large px, incoherent hard parton 
scatterings become dominate and R p a ~ 1 without any 
nuclear modification. The peak in the nuclear modifi- 
cation factor is generally known as the Cronin enhance- 
ment of hadron spectra at intermediate px due to kx- 
broadening through multiple scatterings. Such a peak 
has been observed in pA collisions at the Fermilab ener- 
gies (.Biab = 200 - 800 GeV) |I| and the position of the 
peak at around 4 GeV /c has no apparent energy depen- 
dence and remains the same in HIJING results at both 
RHIC and LHC energies. Nuclear shadowing of parton 
distributions at LHC will only affect the overall value 
of R p a at intermediate px, but the peak- feature and its 
position remain the same. This is a unique prediction 
of the fc^-broadening due to multiple scattering as im- 
plemented in HIJING2.1. It is in sharp contrast to the 
prediction of disappearance of the peak in gluon satura- 
tion model or relocation to much higher px due to strong 
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FIG. 5: (color online) Nuclear modification factor for charged 
hadrons in p + A collisions from HIJING 2.1 with different 
options: default (solid), DHC (dashed), DHC without parton 
shadowing (dot-dashed) and DHC without parton shadowing 
and independent fragmentation (dotted), as compared to data 
from the PHENIX and STAR experiment [H, \2^. The ar- 
rowed lines indicate the most possible behavior of the nuclear 
modification factors at LHC between low and high pr regions. 



anti-shadowing in the collinear pQCD parton model. Ex- 
perimental test of these features will shed light on the 
parton dynamics in high-energy nuclei. 

At very large px, parton shadowing in nuclei, soft-hard 
coupling and enhanced gluonic jets due to valence quark 
conservation all lead to suppression of charged hadron 
spectra. Fragmentation of the projectile string system 
with multiple gluon jets in p + A collisions seems to fur- 
ther suppress high-py hadron spectra as compared to the 
independent fragmentation (dotted). These nuclear ef- 
fects are consistent with existing d + Au data at RHIC. 
However, p + Pb collisions at LHC will be able to test 
these scenarios within the available px range. Since par- 
ton shadowing will disappear at large px [J] due to QCD 
evolution and hard scatterings are de-coherent from soft 
interactions, the nuclear modification factor at LHC will 
likely follow the default and DHC results at low px and 
approach DHC without shadowing at large px with possi- 
ble further modifications due to hadronization of multiple 
jets. In Fig. we indicate this trend by joining these two 
set of results with arrowed lines at intermediate px 

The initial-state interaction and fragmentation of mul- 
tiple jet systems should also affect the final hadron spec- 
tra in high-energy heavy-ion collisions in addition to jet 
quenching due to final-state interaction. The existing ex- 
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FIG. 6: (color online) Nuclear modification factors for final 
parton (upper panel) and charged hadrons (lower panel) in 
central Pb+Pb collisions at ^ = 2.76 TeV from HIJING 2.1 
with different options: default (solid), DHC (dashed), DHC 
without parton shadowing (dot-dashed) and DHC without 
parton shadowing and independent fragmentation (dotted), 
as compared to data from the ALICE and CMS experiment 
[27L I2H ]. Jet quenching due to final-state interaction is not 
considered in the HIJING calculations. 



pcrimental data on hadron spectra modification should 
provide some constraints on the initial-state parton inter- 
actions and the final hadronization mechanism. Shown 
in Fig. [5] are nuclear modification factors for final par- 
ton (upper panel) and hadron spectra in central Pb + Pb 
collisions at yfs = 2.76 TeV. Note that there is no jet 
quenching effect in the calculation shown. On the par- 
ton level, the nuclear modification at large px mainly 
comes from parton shadowing which should be small if 
the scale dependence of shadowing is taken into account. 
The final hadron spectra at large pt are however signifi- 
cantly suppressed due to the Lund string fragmentation 
of multiple parton systems (solid and dashed curves in 
the lower panel) leaving no room for additional suppres- 
sion by jet quenching as compared to experimental data 
from ALICE and CMS @]. This is a strong indica- 
tion that the Lund string fragmentation for hadroniza- 



tion is no longer a valid mode for multiple gluonic jet 
systems in heavy-ion collisions where a deconfined quark- 
gluon plasma is formed. Recent studies show that mul- 
tiple parton scatterings inside a quark-gluon plasma can 
in fact change the the color flow of the jet parton shower 
from that in vacuum (2914321 ] . One therefore should ex- 
pect the fragmentation processes in heavy-ion collisions 
to be significantly modified from that in vacuum. An 
independent jet fragmentation mechanism such as par- 
ton recombination model might be a better alternative 
where local parton phase distributions determine the fi- 
nal hadron spectra and global color flow of the multiple 
jet system is no longer important. Since multiple gluonic 
jet systems in pA collisions can be considered as pre- 
cursors of the quark-gluon plasma, investigation of the 
final hadron spectra at large pt in p + A collisions can 
shed light on the hadronization mechanism of energetic 
jets after they escape jet quenching and emerge from the 
quark-gluon plasma in high-energy heavy-ion collisions. 

In summary, we have studied the nuclear modification 
of hadron spectra in p+A collisions at both the RHIC and 
LHC energies within HIJING2.1 Monte Carlo model and 
their sensitivities to cold nuclear effects. Multiple parton 
correlations in the projectile due to valence quark num- 
ber conservation, string fragmentation of multiple gluon 
jets as well as parton shadowing can all lead to suppres- 
sion of high-pT hadron spectra in p + A relative to p + p 
collisions. Comparisons between HIJING2.1 results with- 
out jet quenching and existing data on large pt hadron 
suppression at LHC indicate an independence fragmen- 
tation mechanism for multiple jet systems such as in pA 
and AA collisions. Experimental test of these effects in 
p + Pb collisions at LHC and its implications for hard 
scatterings in A + A collisions is crucial for the study of 
jet quenching and extraction of jet transport properties 
in the hot quark-gluon plasma. 
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